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Abstract 

The typical AdS/QCD models deal with the large- Ac hmit of QCD, 
as a consequence the meson spectrum consists of the infinite number 
of states that is far from the real situation. Basing on introduction of 
anharmonic corrections to the holographic potential, the corrections 
whose existence has been recently advocated, we construct a class of 
bottom-up holographic models describing arbitrary finite number of 
states in the sector of light mesons. Within the proposed approach, 
the spectrum of masses square has the following properties: It is linear, 

~ n, at not very large n, nonlinear at larger n, with the nonlinear 
corrections being subleading in 1/Ac, has a limiting mass, and the 
number of states is proportional to Ac. The considered holographic 
models reflect thereby the merging of resonances into continuum and 
the breaking of gluon string at sufficiently large quark-antiquark sep- 
aration that causes the linear Regge trajectories to bend down. We 
show that these models provide a correct description for the spectrum 
of excited p-mesons. 
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1 Introduction 



Recently the holographic models of QCD have been successfully applied to 
the description of non-perturbative physics of the strong interactions [1-4]. 
For this reason it is interesting to further investigate to what extent such a 
phenomenological approach (often called AdS/QCD) is able to describe the 
low energy QCD irrespective of whether it is related to the original AdS/CFT 
correspondence [5] or not. 

Maldacena's hypothesis [5] provided a promising way to link the fun- 
damental string theories in the low energy approximation to the strongly 
coupled four-dimensional conformal field theories living on the boundary of 
AdS space [6]. As the AdS/CFT method permits to obtain a theoretical 
control over the strongly coupled gauge theories, the next natural step is to 
extend the idea of holographic correspondence to the physical gauge theo- 
ries, such as QCD. This step was realized, albeit speculatively, through the 
AdS/QCD models (see, e.g., [7] for a brief review). The latter approach, 
however, deals typically with the planar limit of QCD, Nc — > oo, i.e. it 
is rather restricted from the very beginning (see discussions in [8] on that 
point). As a consequence, the meson spectrum in the AdS/QCD models 
consists of an infinite number of infinitely narrow states in a sharp contra- 
diction with what we observe experimentally: There are a few of resonances 
in each channel, they have a finite width and the discrete spectrum gradually 
merges into perturbative continuum. In the program of relating the strings 
to the real hadron physics, the next step is therefore called for — the devel- 
opment of holographic models describing a finite number of hadrons merging 
into continuum. In the present Letter, we will construct a class of bottom-up 
holographic models possessing this property. 

The paper is organized as follows. In Sect. 2, we formulate the problem 
and specify the subject matter of this work. Sect. 3 is devoted to the con- 
struction of a general design for the holographic models describing a finite 
number of mesons in the simplest case. Using the modern experimental data 
we show in Sect. 4 that the found form of spectrum seems to be supported 
experimentally. A simple particular realization of the proposed idea is con- 
sidered in Sect. 5. The results obtained, their relation to other works, and 
the possible consequences are discussed in Sect. 6. Sect. 7 concludes our 
analysis. 
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2 Formulation of the problem 



The problem of primary importance in the non-perturbative QCD is the 
description of spectrum of hght hadrons. In this work, we will be interested 
only in the spectroscopy of light resonances, namely we will consider the 
light mesons because their spectrum is well defined in the large- iVc hmit. 
Experimentally, the most investigated sector is that of the p-mesons [9] , for 
this reason we will restrict ourselves to the p-mesons only since only in this 
case there is enough experimental data to make more or less reliable fits. 

The QCD coupling constant is running, this fact makes impossible 
to follow the principles of AdS/QCD correspondence exactly, nevertheless, 
one can attempt to follow these principles as close as possible. According 
to the averaged fit [9], the value of as at the scale of the p-meson mass is 
q;s(0.776 GeV) ~ 0.71, at the scale of the first radial excitation of the p- 
meson q;5( 1.465 GeV) ~ 0.35, at the scale of the highest observed radial 
excitations (2.300 GeV) 0.28. It is seen that in the energy region where 
the most resonances reside the change of ag is moderate (about 20%), on the 
other hand, the theory is strongly coupled (there are resonances), hence, a 
weakly coupled holographic model dual to QCD in that region has chances 
to exist in some approximate sense. Below that region QCD is strong but 
ctg changes rapidly (see, however, the recent indications on the existence of 
infrared fixed point [10] suggesting ^ < 1), above the resonance region the 
two-flavor QCD is approximately scale invariant but is in the weakly coupled 
regime, in both cases the principles of AdS/QCD correspondence are badly 
violated. Thus, we are going to describe holographically the resonance region 
that is restricted by the infrared. Air, and the ultraviolet, Auv, cutoffs. The 
introduction of the UV-cutoff separating the resonance region from the onset 
of the perturbative continuum is crucial for our further analysis and distin- 
guishes our models from other AdS/QCD models existing in the literature. 
The IR-cutoff is fixed by hand in the hard-wall models [1,2, 11] in order to 
model the confinement, in our case an effective IR-cutoff will be introduced 
in a different way. 

Considering QCD with SU{Nc) gauge group, where Nc is finite and arbi- 
trary, we must answer the following question: How many resonances should 
we take into account? The question is not trivial because, generally speak- 
ing, the cutoff Auv depends on Nc- However, there is a simple argument [12] 
that if the spectrum of masses square is linear, ~ p^n (as expected from 
the phenomenology [13] and the semiclassical string models for mesons [14]) 
the number of resonances should be proportional to Nc- Indeed, in the Breit- 
Wigner parametrization of resonances, the width is r„m„ and if this quantity 
becomes equal to the distance between masses square of the neighboring res- 
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onances, r„m„ ~ /x^, we cannot speak any more about resonances as they 
merge into continuum. On the other hand, according to the semiclassical 
string (flux-tube) models [14] r„ ~ mn/Nc, hence, r„m„ ~ m^^jNc ~ jJ^n/Nc 
and we obtain nmax ~ Nc- This expectation will be an important input in 
what follows. 

3 Construction of model 

In order to advance in holographic description of real hadron physics one 
usually extends the existing successful holographic models. There are two 
complementary ways for such extensions, either by introducing new opera- 
tors/fields following the AdS/CFT principles or by modification of the back- 
ground metric. The former seems to be more consistent but at the same 
time more difficult and somewhat ambiguous at the present stage. The latter 
seems rather ad hoc but simpler, for this reason we will exploit this geometric 
way hoping that subsequently the results obtained will be justified from a 
more fundamental point of view with the help of introduction of appropri- 
ate operators/fields. We note in advance that the first step in justifying our 
approach has been already done in [15]. 

The simplest 5d action describing the spectrum of vector mesons is [2] 

S = j <fx^g (^--Lf^^^fM^v^ , (1) 

where g = \ detgMN\ and Fmn = QmVn — OnVm- Generally speaking, from 
the 5d side we must have a Yang-Mills theory with SU{Nf) gauge group [2] 
but for calculation of the mass spectrum one retains the quadratic in fields 
part only so that the Abelian part of Fmn is enough for our purposes. The 
IR boundary condition is that the action is finite at z = oo. The metric is 
parametrized as 

gMNdx^^dx^ = e^^^''^ {dz^ + ri^^dx^dx"), (2) 

with rjfj^i, = diag(— 1, 1, 1, 1). The equation of motion for action possesses 
a solution for the string modes Vm{x, z) which is supposed to be dual to phys- 
ical states of the gauge theory. Fixing the gauge V^z = 0, the corresponding 
equation reads [3] 

d,{e^d,Vn) + q^e\n = Q, (3) 

where Vn must be normalizable solutions for the 4d-transverse components 
which exist only for discrete values of 4d-momentum = m,^. Perform- 
ing the substitution 

Vn = e-^/2^„, (4) 
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Eq. ([3]) takes the form of a Schrodinger equation, 

- + f/(z)^ = mli,^, (5) 

with the potential 

U{z) = \{j^f + \A". (6) 

Thus the spectrum of the 4d gauge theory and the metric of the dual 5d 
theory are related (up to boundary conditions) by the form of potential (E]). 
The fifth coordinate z is known to correspond to the energy scale Q: Q ^ 1/ z. 
Let us divide the potential U{z) into the UV and IR parts, 

f/(z) = ?7uvW + f/iR(^), (7) 
where U{z) > U\jy{z) and U{z) > Ui^{z). 

2— >0 2— >00 

As follows from the AdS / CFT correspondence, the vector wave functions 
must have the UV asymptotics '>p{z) ^ z"^ [11], this dictates U\]y{z) > 

z-'^. The usual choice is U^jy{z) = 3/(4^2) ^j^^t yields e^^^^) = z'"^. The 
conformal isometry of the metric reflects then the conformal behavior of QCD 
in the ultraviolet. The conformal invariance can be broken by two ways — 
either by introducing a hard cutoff zin [1,2,11] (the resulting spectrum m„ ~ 
n does not agree with the phenomenology in this case) or by introducing 
a nontrivial Uir{z). In order to obtain the desired behavior ~ n one 
has to have ?7ir(-2) ~ z"^ at least at z oo, i.e. a potential of the linear 
oscillator type. This idea was realized in the soft-wall models [3] by means 
of introduction of a dilaton field in action ([T]). 

We want to have a spectrum of oscillator type at small n and a finite 
number of discrete energy levels. This is known to require a certain anhar- 
monicity at large z. The general form of the spectrum given by such anhar- 
monic potentials can be derived in a model-independent way by analyzing 
the anharmonic corrections, 

U^^{z) = u'^z^ + az^ + (8) 

where we have taken into account the 0{z^) term because the first anhar- 
monic correction arising from the 0{z'^) term disappears [16]. For the time 
being we neglect the UV contribution U\]\{z) (to be discussed below). 

Considering the anharmonic corrections as perturbations, Eq. ([5]) leads 
to the spectrum [16] 

= 2tu(n + 1/2) -7(n + 1/2)2 + const, ^ = 0,1,2,... (9) 
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where 

It can be seen straightforwardly that one has the spectrum of oscillator type 
at small n if I7I <C u and a finite number of energy levels if 7 > 0. The 
approximate independence of meson masses on Nc imposes uj = 0{N^). 
In order to have the number of states proportional to Nc we must have 
7 = 0{1/Nc), i.e., the nonlinear corrections to the spectrum are subleading 
in the large- iVc counting. 

Consider now the impact of UV region introducing [/uv(^) = ^f^"^, e ^ 0. 
The boundary condition for the oscillator wave functions, i(j{oo) = 0, is then 
supplied by ^(5) = 0, 5 — 0. As a result the even levels in spectrum (jH]) 
disappear, one has n = 1,3, . . . . Replacing n by 

k = n/2-l/2, fc = 0,1,2,..., (11) 

the leading in Nc contribution is given by 

mfc~4w(fc + 3/4) + const. (12) 

Enlarging e leads to a shift up of energy levels (see, e.g., the exact example for 
e = 3/4 in [3]), finally some of the highest levels can vanish. All this, however, 
does not change the general conclusions about the qualitative behavior of the 
spectrum. 

We assume that the general form of potentials should be such that there 
were no possibility for quantum tunneling to the deep infrared, z ^ 00, which 
we do not know how to interpret physically. Together with the requirement 
of finite number of discrete energy levels this leads to the following general 
form of our potentials: They look like a potential well (with the right wall of 
the form z"^) at relatively small z and gradually transform into a plateau at 

large z, Ui-[i{z) > ^uv' -^^S- 1- The asymptotic behavior of the metric 

follows from Eq. ([U]): 

A{z) ^ ±2Auv2 (13) 

(see Eq. (|2])). In order to have the finite action we have to choose the minus 
sign. The number of discrete energy levels is determined by the value of Auv 
above which (plus an arbitrary constant) one has a continuum and by the 
"width" of the potential which we associate with the IR-cutoff Air. As was 
said above, enlarging the factor in front of U\jy{z) one shifts up the energy 
levels, as a result some of the highly excited states can disappear. The 
quantity Air determines (up to a presumably small constant) an effective 
IR-cutoff ziR in the coordinate space: In the vicinity of zir the potential well 
ceases to be of the form z"^ starting to transform into a plateau. 
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Figure 1: The general form of holographic potentials discussed in the present 
paper. The units are conditional. 

4 Comparison with the experiment 

Let us compare our suggestions with the experiment. In the Section "Light 
Unflavored Mesons", the Particle Data [9] reports five p-mesons: p(770), 
p(1450), p(1700), p(1900), and p(2150), the last two states are omitted from 
the summary table but we include them into our fit as they have been inde- 
pendently seen by several collaborations. Taking the experimental masses of 
these resonances with reported errors we can make the least square fit with 
the curve 

ml = An^ + Bn + C (14) 
and estimate the errors in values of parameters, the result is (in GeV^) 

ml ^ (-0.09 ± 0.02)^2 + (L30 ± 0.08)n + (0.71 ± 0.02). (15) 

In the Section "Other Light Unflavored Mesons", the Particle Data [9] 
reports another two p-mesons: p(2000) and p(2270). The former was seen by 
one collaboration only, for this reason we do not consider it. The latter was 
independently seen by two collaborations, its mass is in a good agreement 
with the fit ([IS]). 
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Since the lowest three resonances are particularly well established, it may 
be reasonable to display also the fit based on these three states only, 

ml ^ (-0.37 ± 0.10)722 + (1.91 ± 0.18)n + 0.06. (16) 

Thus the nonlinear correction to the spectrum is in a qualitative agree- 
ment with our analysis both in sign and in magnitude. Having fit ( JTSl) one 
can obtain the numerical values for the parameters of anharmonic potentials. 



5 A toy-model 

To estimate the typical values of input parameters let us consider the follow- 
ing example of potential having the shape "potential well -|- plateau", 

U{z) = (1 - e-^iR (^-^")) ' + C, (17) 

where C is a constant. This is the Morse potential [17] that has been widely 
used for description of vibrations of nuclei (near zq) in a diatomic molecule. 
In our case, by shifting zq > we can imitate Uuy{z) in the region of appli- 
cability of our models, thus the potential ( fT71) is a good interpolation for the 
class of potentials described above. This potential has an advantage of being 
exactly solvable, the corresponding spectrum is [17] 

ml = 2AuvAiR(n + 1/2) - A^,^{n + 1/2)^ + C, (18) 

0<n<^--. (19) 
- - Air 2 ^ ' 

Exploring fit (fTSll we obtain the values of parameter^: Auv ~ 2.34 GeV, 
Air ^ 0.30 GeV, C ~ 0.04 GeV^. One has then seven resonances with the 
masses (in GeV) 

mn = {0.84, 1.39, 1.72, 1.96, 2.12, 2.24, 2.31}. (20) 

According to our previous discussions, the scaling of cutoffs in Nc is Auv = 
C^(a//Vc), Air = 0(l/v^]Vc). It is interesting to note that to the leading order 
in Nc the mass of the ground state is (neglecting a small C) the geometric 
mean of two cutoffs. 



^As seen the numerical value of C is small and compatible with C — Q within the 
experimental errors. For instance, as long as the mass of p(1900) is not well established, 
the value C = is achieved at TOp(i9oo) ~ 1-93 GeV, with the values of cutoffs being 
affected slightly: Auv ~ 2.32 GeV, Air « 0.31 GeV. 
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Potential ffTTj) can be made a representative from our class of potentials 
if we set zq = and add Uuy{z) = ej but the model cannot be then solved 
analytically. At small enough e and assuming Auv ^ Air we obtain in this 
case that the discrete levels reside in the potential well at 2;uv < z < z\^^ 
where ^uv — V^/Auv, -^m — 1/Air. At the choice e = 3/4, the wave 
functions Vn{z), where n is small enough {i.e. one has nearly linear spectrum) 
approximately coincide with those found in [3]. 

6 Discussions 

The proposed class of AdS/QCD models shares the attractive features of 
both hard- and soft-wall holographic models: Like in the hard-wall models, 
there is an (effective) IR-cutoff related to confinement and there is no artificial 
dilaton the physical origin of which is unclear, on the other hand, similarly to 
the soft-wall models, one has approximately linear spectrum and the absence 
of ambiguity in the choice of IR boundary conditions. However, if we include 
the higher-spin fields in the way proposed in [3], the dilaton field seems to 
be inescapabl^ if one wants to have the slope of trajectories independent of 
spin S and the relation ^ ~ n + S" that fits well the known experimental 
data [18, 19]. The shape of the dilaton will be different in our case, up to a 
factor in the exponent it will look like instead of obtained in [3]. 

It should be noted that the coupling constant entering action ([T]) can be 
obtained by matching to the high-energy asymptotics of QCD two-point cor- 
relators in the same way as in [2] (in the case of finite number of resonances, 
the finite energy sum rules [20] should be used). We believe, however, that 
such a matching is performed out the region of applicability of the models 
under consideration. In any case, it is not needed for the derivation of the 
mass spectrum we are concerned. 

An important result of the considered models is that the spectrum con- 
denses at high energies, i.e. the radial Regge trajectories bend down and the 
discrete spectrum ends near the point of zero slope. Within the fiux-tube 
models of mesons, this effect is usually interpreted as the string breaking, 
hence, the proposed models are able to describe this effect: Up to model- 
dependent constants, the string tension in our scheme is (see Eq. ( IT^ ) 

a~ AuvAiR(l-tn-t) = O(Ar0), (21) 

where 

t = AiR/(2Auv) = 0(l/Arj, (22) 

■^Assuming that there are not additional fields in the action ((T]) which are dual to certain 
QCD operators. 
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i.e. 0" is a constant in the large- iVg limit only. 

The issue of non-linear corrections to the string-like spectrum was dis- 
cussed some time ago in the context of QCD sum rules in the planar limit 
(see, e.g., [14,21,22]), in particular such corrections were advocated to be 
exponentially small, e~*" with 6 > 0, in Ref. [21]. If one assumes the scaling 
b = 0{1/ \/W^) and expands the exponential regarding 6 as a small parame- 
ter, the form of exponential correction will be compatible with our results. 

In the minimal version, the considered class of models has 3 parameters: 
The factors in front of the harmonic and anharmonic terms in the holographic 
potential (we can set a = or = in Eq. (IHl) without change of ensuing 
conclusions). This is not very constraining as long as experimentally only a 
few of states are known in each channel. We note, however, that the first 
two parameters are expected to be universal for all channels, hence, adding 
new channels in the model one introduces only one new parameter for each 
tower of resonances — a new arbitrary constant — or even less due to some 
degeneracies [18, 19], thus such models will be much more predictive if many 
kinds of mesons are described simultaneously. 

Recently the necessity of anharmonic corrections to the holographic po- 
tentials was emphasized in Ref. [15] because they reflect holographically the 
contribution of the QCD operators from the Operator Product Expansion for 
the correlators of quark currents [23]. It is interesting to make the following 
observation: The study of the toy-model in Sect. 5 suggests an intriguing 
possibility that the infinite series of anharmonic corrections found in [15] 
might be summed up into some simple function on the 5d side. 

The introduction of UV-cutoff Auv may solve the following problem in 
building the AdS/QCD models: The 5d action is assumed to be local, i.e. 
the higher- derivative terms are supposed to be suppressed. However, there 
is no understanding what scale should suppress those terms. Within the 
AdS/QCD models put forward in the present work, one can speculate that 
the suppression is provided by powers of l/Auy. 

From the physical sense and numerical estimates for the IR-cutoff Air it is 
tempting to associate this quantity with Aqcd as in the hard- wall models [1] . 
According to a usual belief, however, Aqcd is nearly independent of Nc and 
this provides an approximate m„ = 0{N^) scaling for meson masses. In our 
case, the situation is more tricky: Although Air = 0(l/-\/iVc), the scaling 

~ AirAuv = 0{N^) holds due to the existence of the UV-cutoff Auv = 
0{\/Nc), thus we should assume 

^QCD ~ AirAuv. (23) 

If the considered models indeed reflect the real QCD, for the first time we 
are dealing with a situation when we can learn something qualitatively new 
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about QCD from the AdS/QCD models. Specifically, our previous analysis 
suggests the following mechanism for generating Aqcd in the probable ana- 
lytical solution for QCD: One introduces the IR and UV cutoffs and takes the 
limit Air 0, Auv ^ such that the product AjrAuv remains finite and 
should be associated with Aq^^^, that determines the mass scale of the theory 
with massless quarks. Dealing with an effective theory of QCD one works 
with finite Auv ^-nd nonzero Air that are input parameters of the effective 
theory. Precisely this ideology hes behind the models considered in the given 
work. 

7 Conclusions 

In the present paper, we have designed a class of AdS/QCD models that 
describes a finite number of meson resonances whose spectrum is approxi- 
mately Regge-like. Assuming that the number of distinguishable resonances 
depends linearly on N(. we obtained that the non-linear corrections to the 
Regge-like spectrum are subleading in l/Nc- Within the given models, the 
non-linear corrections stem from the anharmonic contributions to the holo- 
graphic potential that determines the mass spectrum. According to the re- 
cent results [15] the anharmonic contributions on the 5d side appear after 
taking into account the QCD operators responsible for the masses of hadrons, 
therefore the design of the considered models is somewhat justified from a 
more fundamental standpoint. The presented approach and that of [15] shear 
also the following important property: The slope of radial trajectories is, up 
to a dimensionless factor of the order of unity, a product of IR and UV cut- 
offs of the model, i.e. the discrete mass spectrum is equally determined by 
both IR and UV sectors of the theory. This interesting feature suggests that 
something similar may happen in the exact solution for the mass spectrum 
of real QCD. 

Acknowledgments 

The work is supported by the Alexander von Humboldt Foundation. I am 
grateful for the warm hospitality by Prof. Maxim Polyakov extended to me 
at the Bochum University. 



11 



References 

G. F. de Teramond and S. J. Brodsky, Phys. Rev. Lett. 94, 201601 
(2005). 

J. Erlich, E. Katz, D. T. Son and M. A. Stephanov, Phys. Rev. Lett. 
95, 261602 (2005). 

A. Karch, E. Katz, D. T. Son and M. A. Stephanov, Phys. Rev. D 74, 
015005 (2006). 

L. Da Rold and A. Pomarol, NucL Phys. B 721, 79 (2005). 

J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. 
Phys. 38, 1113 (1999)]. 

E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998); S. S. Gubser, 
L R. Klebanov, and A. M. Polyakov, Phys. Lett. B 428, 105 (1998). 

H. Forkel, M. Beyer and T. Frederico, JHEP 0707, 077 (2007); 

F. Jugeau, [arXiv:0902.3864 [hep-ph]. 

T. D. Cohen, | arXiv :0805.48T3l [hep-ph]. 

G. Amsler et ai, Phys. Lett. B 667, 1 (2008). 



A. Deur. iarXiv:0901. 21901 [hep-ph]. 

J. Polchinski and M. J. Strassler, Phys. Rev. Lett. 88, 031601 (2002). 

A. A. Andrianov and D. Espriu, Phys. Lett. B 671, 275 (2009). 

A. V. Anisovich, V. V. Anisovich and A. V. Sarantsev, Phys. Rev. D 62, 
051502(R) (2000); D. V. Bugg, Phys. Rept. 397, 257 (2004); E. Klempt 
and A. Zaitsev, Phys. Rept. 454, 1 (2007). 

M. Shifman, in Quark-Hadron Duality and the Transition to pQCD, ed. 
by A. Fantoni, S. Liuti, and O. Rondon (World Scientific, Singapore, 
2006) [hep-ph/0507246] . 

S. S. Afonin. [i[?Xrv:0902.3959l [hep-ph]. 

L. D. Landau and E. M. Lifshitz, Gourse of Theoretical Physics Vol. 3 
(Pergamon, New York, 1965). 

[17] P. M. Morse, Phys. Rev. 34, 57 (1929). 



12 



[18] S. S. Afonin, Eur. Phys. J. A 29, 327 (2006); Mod. Phys. Lett. A 22, 
1359 (2007); Phys. Rev. C 76, 015202 (2007); Int. J. Mod. Phys. A 22, 
4537 (2007); 23, 4205 (2008). 

[19] M. Shifman and A. Vainshtein, Phys. Rev. D 77, 034002 (2008). 

[20] E. de Rafael, *Les Houches 1171 (1997) [hep-ph/9802448] . 

[21] S. S. Afonin, A. A. Andrianov, V. A. Andrianov, and D. Espriu, 
JHEP 04, 039 (2004). 

[22] S. S. Afonin, Phys. Lett. B 576, 122 (2003); Nucl. Phys. B 779, 13 
(2007); S. S. Afonin, A. A. Andrianov, V. A. Andrianov, and D. Espriu, 
Int. J. Mod. Phys. A 21, 885 (2006); Nucl. Phys. B (Proc. Suppl.) 164, 
296 (2007); AIP Conf. Proc. 1030, 177 (2008). 

[23] M. A. Shifman, A. I. Vainstein, and V. I. Zakharov, Nucl. Phys. B 147, 
385, 448 (1979). 



13 



